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^ '■ ABSTRACT 

We examine the spatial and temporal stability of the Hubble Space Telescope's Advanced 
I Camera for Surveys (ACS) Wide Field Camera (WFC) point spread function (PSF) using the two 

square degree COSMOS survey. This is particularly important for studies of weak gravitational 
lensing, where the ability to deconvolve the PSF from galaxy shapes is of paramount importance. 

■ We show that stochastic aliasing of the PSF necessarily occurs during 'drizzling', the image 
I reduction step where individual exposures are combined and the large geometric distortion in the 

■ ACS WFC is removed. This aliasing is maximal if the output pixel scale is equal to the input 
pixel scale of 0.05". Using simulated stars, we show that this source of PSF variation can be 
significantly reduced by choosing a Gaussian drizzle kernel with a size of 0.8 input pixels and 

Qh' by setting the output pixel size to 0.03". This result has been tested mainly on the COSMOS 

Q . images which have four input exposures each dithered by about 60 pixels in the y direction. 

^ I However, the recommendation holds for any ACS WFC data set taken with a similar number of 

. input exposures and a dither pattern that does not cause output pixels to be derived from input 

pixels that are a significant fraction of the ACS WFC chips away from each other in different 
exposures. We show that the PSF is temporally unstable, most likely due to thermal fluctuations 
in the telescope's focus. We find that the primary manifestation of this thermal drift in COSMOS 
images is an overall slow periodic focus change. 
I Using a modified version of the TinyTim PSF modeling software, we create undistorted stars 

in a 30 X 30 grid across the ACS WFC CCDs. These PSF models are created for telescope 
focus values in the range — 10/im to +5/im at one micron increments, thus spanning the allowed 
range of telescope focus values. We then use the approximately ten well measured stars in each 
COSMOS field to pick the best-fit focus value for each field. The TinyTim model stars are then 
used to perform PSF corrections for weak lensing allowing systematics due to incorrectly modeled 
PSFs to be greatly reduced by using this TinyTim method as compared to using stars within the 
COSMOS fields alone, as is typically done in weak lensing analyses. We have made the software 
for PSF modeling using our modified version of TinyTim available to the astronomical community. 
We show the effects of Charge Transfer Efficiency (CTE) degradation, which distorts the object 
in the readout direction, mimicking a weak lensing signal. We derive a parametric correction for 
the effect of CTE on the shapes of objects in the COSMOS field as a function of observation 
date, position within the ACS WFC field, and object flux. Finally, we discuss future plans to 
improve the CTE correction. 

Subject headings: instrumentation: detectors — surveys — techniques: image processing 
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1. Introduction 

The addition of the Advanced Camera for Sur- 
veys (ACS; Pavlovsky et al.2006) to the Hubble 
Space Telescope (HST) in February 2002 enabled 
significant science returns in a variety of disciplines 
due to the increased resolution, areal coverage, 
and quantum efficiency of ACS Wide Field Cam- 
era (WFC) as compared to previous HST imaging 
instruments including WFPC, WFPC2, and STIS. 
Despite the myriad successes of ACS, some areas 
of study have been hampered by a lack of under- 
standing of the properties of the ACS WFC Point 
Spread Function (PSF). In particular, studies of 
weak gravitational lensing, in which the shapes of 
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background galaxies undergo a small coherent dis- 
tortion by foreground dark matter, are made dif- 
ficult by time variability of the ACS PSF. Since 
we are interested in weak gravitational lensing, we 
concentrate our studies on the ACS WFC and do 
not study the PSF of other ACS channels. Here- 
after, when wc refer to the ACS, we arc referring 
to the ACS WFC. Because the slight shape dis- 
tortions due to weak lensing are up to an order of 
magnitude smaller than the shape distortions of 
small galaxies due to the ACS PSF, it is critical to 
have accurate PSF models for use in deconvolving 
the galaxy shapes from the PSF. There arc notable 
exceptions in which ACS has had great success in 
weak lensing, including the study of weak lensing 
by galaxy clusters (e.g Lombardi et al. 2005; Jce 
et al. 2005; Jee et al. 2006), and the results of 
Heymans et al. (2005) using the GEMS survey. In 
the case of lensing by galaxy clusters, the weak 
lensing signal is typically large compared to the 
signal arising from large-scale structure alone and 
thus PSF effects are less important. For clusters 
at z > 1, PSF effects become important, and Jee 
et al.(2006) have found that a limited number of 
PSF models derived from globular cluster observa- 
tions can be used for PSF correction. The GEMS 
survey contained many images taken in a rela- 
tively short amount of time such that the PSF 
time variations could be internally calibrated us- 
ing stars in the survey. However, for surveys taken 
over an extended period of time in which there are 
few stars and the weak lensing signal is low, the 
time variability of the ACS PSF precludes using 
the stars in each exposure to make a PSF model 
suitable for weak lensing PSF deconvolution. In 
this paper we offer an alternative solution in which 
templates of model PSFs are created and the few 
available stars in each ACS field are used to se- 
lect the appropriate PSF template to be used for 
PSF deconvolution. A similar method using tem- 
plates created by combining the stars from many 
exposures with similar PSF properties was used by 
Schrabback ct al. (2006). Another alternative ap- 
proach is to use globualr cluster fields with many 
stars to model the PSF (e.g. Jee et al.2005; Jee 
et al.2006). In addition to the temporally varying 
telescope PSF, degradation of the charge trans- 
fer efficiency (CTE) of the ACS CCDs is gradu- 
ally adding additional variation to the measured 
PSF. The GEMS survey, taken relatively early in 
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the life of ACS did not necessitate a correction for 
this effect but more recent surveys do need to take 
this effect into account. This paper also addresses 
that problem by providing a prescription by which 
measured galaxy shapes can be corrected for the 
effects of CTE degradation. 

The Cosmic Evolution Survey (COSMOS) 2 
square degree field (~ 1.64 square degrees with 
ACS; Scoville et al. 2006) has a unique combina- 
tion of area, depth, and resolution that opens up 
new areas of study with the HST. This is partic- 
ularly true in weak lensing where the surface den- 
sity of resolved galaxies (up to ^ 80 per square ar- 
cminute) is significantly higher than that available 
from ground-based surveys (10 to 30 to per square 
arcminute). Along with the excellent photometric 
redshifts available in the COSMOS field due to ex- 
tensive ground-based follow-up, this allows stud- 
ies of cosmic shear tomography at small angular 
scales to unprecedented accuracy (Massey et al. 
2006a). COSMOS also uniquely affords the oppor- 
tunity for high -resolution dark matter maps over 
a wide field (Massey et al. 2007). Such analyses 
depend critically on the ability to deconvolve the 
effects of PSF from the measured galaxy shapes 
and we demonstrate in this paper a procedure for 
doing so. The demanding nature of our weak lens- 
ing analysis will provide PSF models of immediate 
applicability to other studies, including an exami- 
nation of AGN in the COSMOS field (Gabor et al. 
2006). 

Typical ground-based weak lensing surveys use 
the stars within each field to make a model of 
the PSF. This PSF model, often represented as 
a polynomial across the field, is then used to de- 
convolve the PSF from galaxy shapes using a va- 
riety of methods (see Heymans et al. 2006 and 
Massey et al. 2006b for a description of many of 
these methods). The modeling and deconvolu- 
tion of the PSF remains the primary systematic 
for weak lensing studies and considerable work is 
currently going into optimizing the methods for 
doing this (again, see Heymans et al.2006). HST 
images useful for weak lensing contain too few 
stars in each image to make a viable PSF model. 
In the past, this has been overcome by collecting 
stars from many images to model the PSF (e.g 
Rhodes, Refregier, & Groth 2001; Hoekstra et al. 
1998). However, we find that the ACS PSF is 
not sufficiently stable to allow for the collection of 



stars from many images that will provide models 
at the accuracy level we needed for cosmic shear. 
Thus, we have modified the PSF creation software 
TinyTim (Krist & Hook 2004) to create simulated 
PSFs that we use for PSF correction. The rea- 
sons the ACS PSF is more temporally unstable 
than that of previous cameras such as WFPC2 
and are two-fold. First, WFPC2 is on-axis , while 
ACS is significantly off-axis yielding a more sensi- 
tive measure of defocus. Second, ACS has a more 
asymmetric distortion compared to the largely ra- 
dially symmetric pattern for WFPC2. These ef- 
fects combine to provide more sensitivity to focus 
changes. 

The degradation of CTE is caused by charge 
traps created by charged particle radiation dam- 
age induced CCD defects (Riess & Mack 2004; 
Mutchler & Sirianni 2005). These traps cause 
a trailing of the charge during readout that can 
mimic a weak lensing signal. This effect is partic- 
ularly insidious for weak lensing applications be- 
cause it preferentially affects low flux objects (such 
as galaxies) rather than high flux objects (such as 
bright stars). So, the typical weak lensing method 
of using high signal-to- noise (S/N) stars to correct 
the PSF of low signal to noise galaxies completely 
fails in the case of CTE degradation. Nor, strictly 
speaking, does it act as a convolution, although 
it has been treated as such to first order in the 
past (e.g. Rhodes et al. 2004). The effect of CTE 
degradation depends in a complicated way on the 
date of observation, the position within the CCD, 
the fiux of a source, and the background level of 
the image. The optimal solution to the problems 
caused by CTE degradation would be to under- 
stand and remove the effects of CTE on the im- 
ages themselves as the first step in image reduc- 
tion. However, the effects of CTE on individual 
pixel values has not yet been fully quantified and 
will be investigated in a future paper. We cor- 
rect for CTE effects by subtracting a model of the 
spurious (CTE induced) signal from the measured 
galaxy ellipticities at the end of our weak lensing 
pipeline. 

This paper is organized as follows. §[2] describes 
the data and software tools we use to analyze the 
data. This section includes descriptions of the 
COSMOS survey and auxiliary data sets we use 
in § 12.11 a brief description of the multidrizzlc im- 
age reduction pipeline in § 12. 2| a description of the 
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TinyTim PSF simulation software in § 12.31 and 
a description of our weak lensing PSF correction 
scheme in ^2.4\ We detail the image reduction pro- 
cedures we use for weak lensing images in § [3) In 
§ m we describe the Tiny Tim-based PSF models 
we have created and the procedure we use to de- 
convolve the PSF from galaxy shapes is given in 
§[3 We conclude in §[6l where we describe how the 
methods we have detailed can be of use on other 
data sets and with other weak lensing methods. 
In the conclusion we describe steps we have taken 
to make our PSF modelling software available to 
the community for use on other HST data sets. 

2. Data and Software 
2.1. Data 

The primary motivation for our analysis of the 
ACS PSF is to perform various weak gravitational 
lensing analyses on the COSMOS field. The COS- 
MOS field is a contiguous square covering 1.637 
square degrees centered at 10:00:28.6, +02:12:21.0 
(J2000) COSMOS was imaged with the ACS Wide 
Field Camera (WFC) during during HST cycles 12 
(proposal ID 9822) and 13 (proposal ID 10092), 
between October 2003 and June 2005. This data 
set contains 575 slightly overlapping pointings, 
each with four individual exposures of 507 sec- 
onds dithered by (5,60) pixels. At this depth, each 
10 square arcminute COSMOS pointing contains 
approximately 800 galaxies useful for weak lens- 
ing (sufficient size and signal-to-noise (S/N) and 
approximately 10 stars suitable for measuring the 
PSF. For a more complete description of the COS- 
MOS field see ScoviUe et al.2006. COSMOS was 
imaged in the F814W (/) filter and our measure- 
ments of PSF are specific to that filter. However, 
the methods we have developed are general and 
can be applied to other ACS WFC filters using 
the software we have made public (see S}6]). The 
applicability of our software and methods is lim- 
ited by the ability of our modelling software, which 
does not include such effects as the long wave- 
length (e.g. red) scattering within the ACS WFC 
CCDs (see for instance, Sirianni ct al.2005). 

The ACS WFC's PSF varies due to focus 
changes caused by thermal fluctuations of the 
HST on orbit (Heymans et al. 2005; Rhodes et al. 
2006; Schrabback et al. 2006; Lallo et al. 2005; 
Lallo et al. 2006; Makidon et al. 2006a; Makidon 



et al. 2006b). This manifests itself in PSF pat- 
terns that vary significantly over time (Figure [T]). 
Since the changes in focus are stochastic, the focus 
of a particular observation cannot be predicted a 
priori and must be inferred from the data. We 
have developed a method for determining the tele- 
scope focus (see ^ and we show in Figure [2] the 
focus of HST as a function of time over the course 
of the COSMOS observations. There are two time 
scales over which the HST focus changes. The 
first time scale is the 90 minute HST orbit; the 
telescope expands and contracts, or "breathes," 
as the telescope goes into and out of sunlight (see, 
for instance, Lallo et al.2006). Since each COS- 
MOS observation takes place over a single orbit, 
we average over this source of focus variation in 
COSMOS (and indeed in most HST apphcations). 
This intra-orbit breathing is a large source of the 
individual error bars on the points in Figure [2] be- 
cause the intra-orbit focus changes are of order a 
few microns or more (Lallo et al. 2006). There is 
also an apparent slow oscillatory drift in the focus 
of the telescope over several weeks; it is this grad- 
ual drift that we attempt to model and correct 
for as described below. From Figure [2] it is clear 
that sampling the focus every month as done by 
Lallo et al. provides insufficient sampling to satis- 
factorily model this drift. We find that the focus 
is typically within ±2^m of — 3/im, as shown in 
Figure [H During the first decade of HST oper- 
ations it was necessary to periodically (approxi- 
mately every 6 months) reset the telescope focus 
as the struts holding the secondary mirror in place 
slowly outgassed water and shrank. However, the 
last commanded focus changes of HST were in 
December 2002, December 2004, and July 2006 
because this slow focus drift due to the shrinking 
of the structure has stabilized. 

In order to examine the behavior of the ACS 
PSF we make use of two other HST data sets. The 
first are images of the stellar field 47Tuc taken as 
part of the ACS calibration program (GO-10048; 
Krist 2003). The three 47Tuc images are single 
30 second exposures of dense stellar fields taken 
between 28 November 2003 and 7 September 2004 
with the F814W filter. These images are processed 
with the same pipeline we describe in § [3] and we 
use these images to test our PSF models and focus 
prediction methods as described in § |4l 

We also use Extended Pixel Edge Response 
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Fig. 1. — These figures demonstrate that the 
ACS PSF is time variable. Each tick mark repre- 
sents the magnitude and direction of the eUiptic- 
ity of a star in the COSMOS field. The left panel 
shows stars for the HST near the "nominal" focus 
value and the right panel shows stars in images 
taken when the primary-secondary mirror separa- 
tion was 3 microns closer than nominal. Since the 
focus varies with time, the PSF is not temporally 
stable. While the plots are noisy, the different un- 
derlying patterns can be clearly discerned. 
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Fig. 3. — A histogram of the best-fit focus values 
for each of the 575 COSMOS fields. 
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Fig. 2. — The top panel shows the focus of HST 
for each field in the COSMOS survey(small data 
points) . The errors bars are the error on the mean 
focus determined from the ~ 10 useful stars in 
each COSMOS image. The larger, fainter data 
points are focus values determined from a moni- 
toring program using the ACS HRC (Lallo et al. 
2006). Their error bars represent the full focus 
variation during an HST orbit. The date of the 
2004 secondary mirror movement is labelled. The 
bottom panel shows the focus for a small time win- 
dow during Cycle 12 (March to May 2004). The 
cyclical pattern of the focus is clear in the bottom 
panel. We describe how we determine the tele- 
scope focus value in fj4] 



(EPER) data from HST calibration program 
10369, consisting of images acquired using inter- 
nal targets (lamps) only, taken during Earth oc- 
cultation time (Mutchler & Sirianni 2005). These 
images overscan the ACS CCD fiducial range of 2 
X 2048 X 4096 pixels, providing additional 'vir- 
tual pixels' along the edges of the chips, into which 
pixel charge (but not actual signal from the lamps) 
can transfer during readout. The signal that oc- 
curs in the overscan region is charge that is re- 
leased from charge traps during the readout pro- 
cess. This allows calibration of the CTE of the 
ACS CCDs. The EPER images are 1.2 second ex- 
posures, taken using a combination of the F555W 
and F435W filters to reduce the incident power 
from the lamps. Note that the CTE occurs within 
the CCDs so is completely independent of the fil- 
ters and only depends on flux, date, and CCD 
position (see i j5.2[) . 

2.2. Multidrizzle 

The ACS is situated off-axis in the HST fo- 
cal plane and the ACS focal plane is not nor- 
mal to incident light rays. Therefore ACS im- 
ages have a large geometric distortion (Anderson 
2006; Pavlovskyet al.2006). The 4096 x 4096 ACS 
WFC CCDs do not sample a square on the sky. 
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they sample a "squashed parahelogram" (see, for 
example, Figure [H The program MultiDrizzle 
removes this geometric distortion while simulta- 
neously removing cosmic rays and bad pixels as 
well as combining multiple exposures into a sin- 
gle output image (Koekemoer et al. 2002). The 
geometric distortion is easily measured and cor- 
rected for with astrometric observations of stellar 
fields, as shown by Meurer et al. (2002) and An- 
derson (2006) . They fit the geometric distortion to 
a quartic polynomial and the resulting fit is good 
to better than 0.05 pixels across the ACS WFC. 
Subsequent observations have improved this fit 
and the latest distortion polynomial coefficients 
are readily available from the Space Telescope Sci- 
ence Institute. Despite the fact that the geometric 
distortion is the largest PSF effect, it is the easi- 
est to correct for and we use the default correction 
available via MultiDrizzle. The geometric dis- 
tortion is not completely temporally stable and, 
indeed, has changed over the lifetime of the ACS 
(Anderson 2006). However, the size of this varia- 
tion is much smaller than the variation in PSF due 
to changes in the telescope's focus. The latest geo- 
metric distortion corrections are used in the reduc- 
tion of the COSMOS data (Koekemoer et al.2006). 
The pixel scale of the MultiDrizzle output image 
can be smaller than the pixel scale of the input im- 
age (or images). The details of how MultiDrizzle 
is run on COSMOS data are given by Koekemoer 
et al. (2006). We make some important changes 
in the default MultiDrizzle parameters in order 
to optimize the images for weak lensing analysis 
as described in § [3l 

2.3. TinyTim 

There are not enough stars in each extragalac- 
tic COSMOS ACS image to allow us to model the 
PSF across the field through interpolation. There- 
fore, we have simulated artificial stars in ACS 
fields at arbitrary positions (in the following dis- 
cussion we will use the term 'star' and 'PSF' inter- 
changeably) using the TinyTim software package 
(Krist & Hook 2004). TinyTim can create PSFs for 
any current HST camera with any filter combina- 
tion, at any detector position, and for any given in- 
put spectrum. TinyTim creates FITS images con- 
taining stars that include the effects of diffraction, 
geometric distortion, and charge diffusion within 
the CCDs. By default, the stars appear as they 



would in raw HST images. In the case of ACS, 
this means the stars are highly distorted and have 
a pixel scale of 0.05 arcseconds per pixel. TinyTim 
is able to make highly oversampled PSFs and can 
incorporate the focus position (primary/secondary 
mirror spacing of the telescope) . 

We have adapted version 6.3 of the TinyTim 
software package to create simulated ACS starfields. 
By default, ACS creates only single stars or several 
stars in a small portion of the ACS WFC. Using 
the IDL programming language, we have written 
a wrapper that allows us to run TinyTim multi- 
ple times and create a grid of PSF models across 
the whole ACS field of view. This is available from 
\protecttittp://www. astro. caltech.edu/~rjm/acs. 
We insert our artificial stars into blank images 
with the same dimensions and FITS structure as 
real ACS data, thereby manufacturing arbitrar- 
ily dense starfields. The default TinyTim pipeline 
calculates a diffraction pattern (spot diagram) , ge- 
ometrically distorts this pattern, and adds charge 
diffusion between adjacent pixels. The geometric 
distortion coefficients built into TinyTim were cur- 
rent at the time of the program's release but have 
since been superseded. Our IDL TinyTim wrapper 
allows us to input the most recent distortion coef- 
ficients. These three effects occur (in that order) 
in real data and depend upon the position of the 
star in the ACS field of view. TinyTim is run in 
steps, and, with our modifications, the TinyTim 
software allows some of these steps to be skipped. 
By skipping various steps within our IDL wrapper 
we allow for two important changes to this basic 
pipeline. The deviations from the default TinyTim 
pipeline are: 

• In order to examine the effects of the dis- 
tortion removal process (MultiDrizzle in 
our case), we incorporate the option to al- 
low each star to have an identical diffraction 
pattern and charge diffusion, but a geomet- 
ric distortion determined by the location of 
the PSF within the ACS field of view. Once 
the geometric distortion is removed by run- 
ning the field through MultiDrizzle these 
stars should all appear identical. 

• We allow TinyTim to create starfields that 
do not contain the effects of geometric dis- 
tortion at all, instead modelling stars as they 
would appear after a perfect removal of geo- 
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Fig. 4. — A typical artificial starfield created using 
our modified version of TinyTim. This field is at 
focus of — 2/ini. 

metric distortion. Conversion between dis- 
torted and non-distorted frames, which is 
necessary to simulate charge diffusion in the 
raw CCD, is performed using very highly 
oversampled images. This avoids stochastic 
aliasing of the PSF (see Q, and minimizes 
noise in the PSF models. This makes use 
of the geometric distortion coefficients de- 
scribed in Gonzaga et al. (2005). 

Our modified version of TinyTim thus allows 
us to make dense artificial starfields at a range of 
focus values, and is publicly available from 
pittp://www. astro. caltech.edu/~rjm/acs We have 
created starfields in the F814W filter for the COS- 
MOS analysis, but the procedures we have devel- 
oped are general and can be adapted to any filter. 
For reasons of speed, we raytrace through the op- 
tical path at a single wavelength of 814nm. Using 
a full stellar (or galaxy) spectrum in the F814W 
filter does not significantly change the output PSF 
but adds considerably to the processing time. We 
generate the PSF models in an oversampled pixel 
frame. This can be repixelated to the ACS pixel 
size in geometrically distorted (raw) coordinates 



or in geometrically corrected coordinates. In the 
distorted coordinates, we find that the simulated 
lightly smaller than real stars, possibly 
because of pixelization effects. This is fixed by 
convolving the PSF model with a square kernel 
the same size as the pixels. Note that, as dis- 
cussed in ^2.1\ our TinyTimPSF models do not 
include all PSF effects, especially near the red end 
of the spectrum. 

We also have included the ability to create a 
set of raw images with an arbitrary dither pat- 
tern. Stars can be added to each geometrically 
distorted dithered image in such a way that all 
the stars line up with each other on the output im- 
age. This portion of our TinyTim pipeline makes 
use of the "wtranback" coordinate transformation 
routine built into pyraf . We can then make four 
dithered images with, for instance, the COSMOS 
dither pattern and then combine them into a single 
output file using MultiDrizzle. See, for example, 
Figure [5] This is non-trivial because the large geo- 
metric distortion in ACS can substantially change 
the relative spacing of stars in dithered images. 
We describe the results of such tests in the follow- 
ing section. 

2.4. Shape Measurement 

Here we introduce our formalism for the mea- 
surement of galaxy shapes. The method of 
Rhodes, Refregier, & Groth (2000; hereafter 
RRG) has been optimized for space-based images 
with small PSFs and has previously been used 
on weak lensing analyses of WFPC2 and STIS 
data (Rhodes, Refregier, & Groth 2001; Refregier, 
Rhodes, & Groth 2002; Rhodes et al. 2004). Fol- 
lowing the formalism of RRG, we parameterize ob- 
ject shapes by measuring the Gaussian weighted 
second order moments: 
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The summation is over all pixels, w is the size of 
the Gaussian weight function, / is the pixel inten- 
sity, and the coordinates measured in pixels. 
These moments are used to derive the ellipticity 
of an object Ci and a size measure d given by 
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The RRG method also requires the measurement 
of the five fourth order moments given by: 



(5) 



in order to correct for the efi^ects of the weight- 
ing function and the shape distortions caused by 
the PSF. These fourth order moments are used to 
make small corrections to the second order mo- 
ments. In this way, the RRG method is a bridge 
between the earlier KSB (Kaiser, Squires and 
Broadhurst 1995) method which uses only second 
order moments, and more advanced methods such 
as Shapelets (Massey & Refregier 2005) which uses 
moments to arbitrarily high order based on the 
amount of information available in each object. 

The ellipticities are related to the shear, the 
quantity of interest for weak lensing, via the shear 
susceptibility factor G: 



(6) 



Previous weak lensing analyses using the RRG 
method have made use of a single value of G for the 
entire survey. We find that residual PSF system- 
atics are reduced when G is allowed to vary as a 
function of object size d, ellipticity e and signal-to- 
noise. Leauthaud et al. (2006) contains a discus- 
sion of the calculation of G for the COSMOS data 
set. Given an output pixel size of 0.03 arcseconds 
(see Sj3]and the Source Extractor parameters in 
Leauthaud et al.(2006), we find that the optimal 
weight function size for COSMOS data is 



w 



maximum[2v ab, 6] 



(7) 



where a and b are the Source Extractor com- 
puted semi-major and semi-minor axes, and the 
minimum weight function width 6 has been empir- 
ically determined to be the optimal weight func- 
tion to measure the shapes of stars. These above 
parameters (especially the factor of 2 in Equa- 
tion [7]) have been tuned empirically and depend 
on the Source Extractor settings. 



3. Optimal Image Reduction 

The transformation of pixels from a distorted 
input image to an undistorted output plane can 
introduce significant "aliasing" of pixels if the out- 
put pixel scale is comparable in size to the in- 
put scale. When transforming a single input im- 
age to the output plane, point sources can be en- 
larged, and their ellipticities changed significantly 
depending upon their sub-pixel position. This is 
one of the fundamental reasons why dithering is 
recommended for observations, since the source is 
at a different sub-pixel position in different expo- 
sures, thus partially mitigating these effects when 
exposures are combined. However, the residual 
aliasing in combined images is sufficient to pre- 
vent the measurement of small, faint galaxies at 
the precision required for weak lensing analysis. 

In order to demonstrate the effects of aliasing 
we created a set of TinyTim starfields that con- 
tain the same diffraction and diffusion for each star 
but a geometric distortion given by the position of 
the star within the ACS WFC field. We then use 
MultiDrizzle (with default settings) to remove 
this geometric distortion. The resulting, undis- 
torted field should contain stars that are identical 
in all portions of the chip. We show in the first 
panel of Figure [5] that this is not the case. In this 
figure, each tick mark represents a TinyTim cre- 
ated star. The length and orientation of the tick 
mark represents the size and direction of the star's 
ellipticity as given by Equations [2] and [3l The fact 
that the tick marks vary is evidence of aliasing. 
We have found that dithering and reducing the 
size of the output pixels reduces the scatter in the 
ellipticities of stars in this study. The largest gain 
comes from reducing the size of the output pix- 
els, but, as expected, the scatter in e also drops as 
roughly ^/N for TV dithers. This confirms the idea 
that the repixelization adds stochastic noise to the 
ellipticity when the sub-pixel positions are uncor- 
related. Thus, sub pixel aliasing is not a problem 
when there are many dithers (e.g. the Hubble Ul- 
tra Deep Field) but for fields with a limited num- 
ber of dithers like COSMOS, this is a significant 
source of PSF error. As we show below, this ef- 
fect can be minimized by carefully choosing several 
MultiDrizzle parameters. 

Pixelization effects are unavoidable during the 
initial exposure, when the detector discretely sam- 
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2000 4000 6000 SOOO 

Single Image, 0,03 arcsecond pixels" 

Fig. 5. — The "aliasing" of the PSF introduced 
when transforming distorted input images to an 
undistorted output frame. Each tick represent 
the elhpticity of a star that has undergone iden- 
tical diffraction in the telescope's optics. Each 
tick mark should therefore look identical. The 
difference between stars is the sub-pixel posi- 
tion when geometric distortion is removed via the 
MultiDrizzle program. The first panel shows 
a single image run through MultiDrizzle with 
the default settings. This problem can be ame- 
liorated by altering several of the MultiDrizzle 
settings and using dithered input images. The 
second panel shows that the ellipticity scatter is 
reduced when four dithered images are combined 
with the default MultiDrizzle settings. We show 
in Figure [7] and in the third panel that reducing 
the output pixel scale greatly reduces this source 
of PSF uncertainty. This panel shows the re- 
sults when a single input image is processed us- 
ing our optimized MultiDrizzle settings, includ- 
ing a final output pixel scale of 0.03 arcesconds. 
The scatter in the output ellipticities is further re- 
duced by using 4 dithered images with optimized 
MultiDrizzle settings. We do not show this plot 
as the reduced ellipticity scatter is difficult to de- 
tect by eye. The inter-chip gap visible in the first 
and third panels is covered by the dither pattern, 
which was one of the goals of the COSMOS dither 
pattern. 



pies an image. However, it is clearly desirable 
to minimize related effects during data reduction. 
The effect on each individual object depends on 
how the input and output pixel grids line up. 
These arguments are demonstrated graphically in 
Figure [6l That figure demonstrates that we can 
mitigate this effect by using a finer grid of output 
pixels {e.g. Lombardi et al. 2005). This reduction 
in pixel scale (which will cause a corresponding in- 
crease in computer overheads) can be performed in 
conjunction with simultaneously "shrinking" the 
area of the input pixels that contains the signal, 
by making use of the MultiDrizzle pixf rac pa- 
rameter and convolution kernel. 

We have run a series of tests on the simu- 
lated PSF grids to determine the optimal values 
of the MultiDrizzle parameters specifically for 
weak lensing science. As described above, we pro- 
duced a grid of stars that ought to look identi- 
cal after the removal of geometric distortion. We 
then ran a series of tests using MultiDrizzle on 
the same input image but with a range of out- 
put pixel scales, convolution kernels, and values 
of pixfrac. We measured the scatter in the el- 
lipticity values in the output images. The smaller 
the scatter, the more accurately the PSF is rep- 
resented. We found the results were not strongly 
dependent on the choice of pixfrac and settled 
on pixfrac=: 0.8. We show in Figure [7] that PSF 
stability is improved dramatically by reducing the 
output pixel scale from 0.05 arcseconds (the de- 
fault) to 0.03 arcseconds. There is a very slight 
gain in going to smaller output pixel sizes but the 
storage requirements increase rapidly as the num- 
ber of output pixels is increased. The gain in go- 
ing to smaller pixel scales is more stable with a 
Gaussian kernel than with the default square ker- 
nel. Therefore, for weak lensing we use an output 
pixel scale of 0.03 arcseconds, pixf rac=0.8, and 
a Gaussian kernel. Despite its clear advantages 
for weak lensing studies, the Gaussian kernel does 
have some general drawbacks, such as the intro- 
duction of more correlated noise which may not 
be desirable for other types of science where min- 
imization of correlated noise is important. 

4. Focus Dependent PSF models 

Changes in the HST's focus significantly change 
the PSF. These focus changes are caused by a 
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A first pixelization at the focal plane is unavoidable 
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A second coarse pixelization merely adds shape noise 
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A second pixelization with higher output resolution, and a kernel 
that better preserves the action of a PSF as a pure convolution 
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Fig. 6. — This series of illustrations depicts the aliasing of images caused by pixelization. The first aliasing 
effect (depicted in the top two rows) is an inevitable consequence of observing the continuous sky with 
discrete pixels. The top row of this panel shows a simulated object on a pixel grid, and the second row shows 
how that object would appear in a pixelized image. Point sources that lie near the center of a pixel are 
mostly detected in that pixel, as shown in the center column. However, if the center of an object lies near 
the border between two pixels, the object is elongated in one direction during the process of observation. 
Such an elongation mimics the shear caused by weak lensing and is shown in the left and right columns 
of the first two rows. The middle tow rows illustrate a second pixelization that can change the shapes of 
observed objects. This pixelization occurs during the MultiDrizzle stage of image processing when images 
are resampled onto a new output grid in order to remove geometric distortion. In each of the three instances 
here, a perfectly centered star in the original (black tilted) coordinate frame is drizzled onto the underlying 
output pixel grid, which is aligned with the page. A second aliasing occurs at this stage if the input pixel 
lies on the border of two output pixels. Again, the middle columnn shows no aliasing but the left and right 
images show the object being elongated in one direction. The image distortion in the second pixelization is, 
however, somewhat avoidable. If the output pixel grid is made smaller, as shown in the bottom two rows, 
the shape of the object on the final pixel grid is much less dependent on the relative alignments of the input 
and output pixels. Sub-pixel dithering of many input images and subsequent image combination onto a fine 
output pixel grid can further reduce the effects of pixelization. 
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Fig. 7. — RMS ellipticity introduced during 
MultiDrizzle. Lower values of the scatter in el- 
lipticity show more stable behavior of the PSF 
during this stage of image combination and ge- 
ometric distortion removal. The only sources of e 
variation here is pixelization and measurement er- 
rors; each star was designed to have the same PSF 
after the removal of geometric distortion. Each 
measurement represents a grid of 20 x 20 stars 
spread across the ACS WFC field. The scatter 
in the ellipticity is measured as the standard de- 
viation of the mean of measured ellipticity val- 
ues. Note the resonance in the square kernel at 
0.025 arcsecond pixels, half the original pixel size. 
The reduction of output pixel scale is the most 
important improvement we make on the default 
MultiDrizzle pipeline, with a slight further gain 
from going to a Gaussian rather than square ker- 
nel. For the COSMOS weak lensing images, we 
use a Gaussian kernel and an output pixel size of 
0.03 arcseconds (Leauthaud et al.2006). This min- 
imizes the effects of undersampling on the PSF 
and produce images that are optimal for weak 
lensing science at the cost of introducing signif- 
icant additional correlated noise relative to the 
square kernel. This is not critical for weak lensing 
science but may not be optimal for other types of 
science. 

change in the primary-secondary mirror spacing 
brought about by thermal fluctuations. This spac- 
ing (which we hereafter refer to as the focus) can 
deviate from its nominal value in the range — 10/im 



to +5nm. We thus want to know what the PSF 
looks like at each point in the ACS WFC field 
for this entire range of focus values. The repix- 
elization of ACS data necessary to remove geo- 
metric distortion causes stochastic aliasing of the 
PSF even with the optimal MultiDrizzle param- 
eters presented in § [3l Because of this, we cre- 
ate TinyTim stellar fields without geometric dis- 
tortion. These contain stars that have diffraction 
and diffusion given by their position within the 
ACS WFC field, but appear as they would after a 
perfect removal of geometric distortion. These are 
stars as they would appear if the number of expo- 
sures N was very large. Since we are creating sim- 
ulated stars, we make the spacing between stars 
small enough that we do not have to do a compli- 
cated interpolation between stars to find the PSF 
at a particular point in the ACS WFC field. We 
simply choose the closest star. We find that this 
condition is satisfied if we create grids of 30 x 30 
stars across the two ACS WFC CCDs. We also 
find that creating such a grid in the focus range 
— 10 to +5fim at one micron increments gives suf- 
ficient resolution in focus that the difference in el- 
lipticity between a star at two adjacent focus val- 
ues (at the same chip position) is smaller than the 
measurement error in measuring the ellipticity of 
a star in a typical noisy COSMOS image. 

In a typical COSMOS field approximately 10 
stars have a suitable S/N and are sufficiently de- 
blended from other objects to provide accurate 
PSF measurements. See Leauthaud et al. 2006 
for a description of star-galaxy separation in the 
COSMOS weak lensing catalog. We compare these 
stars to each focus model in the range [—10, 5]/im. 
We determine the TinyTim focal position that 
minimizes the between the ellipticities of the 
stars in data and the corresponding closest stars 
in the TinyTim model, where the is defined as 

x' = Y.(^i~err + {e;-err («) 

The superscripts * and TT represent real COSMOS 
stars and TinyTim simulated stars, respectively, 
and the sum is over all the stars in the image that 
make the cuts described above. Thus, we deter- 
mine the best fit for the telescope's focus at the 
time of the observation. This is the focus value we 
show in Figure [2l Note that this method of deter- 
mining focus necessarily averages over any intra- 
orbit focus changes. Figure [8] shows the stellar 
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ellipticity pattern in the ACS WFC for a value of 
focus=— alongside the pattern formed by av- 
eraging all of the approximately 2000 COSMOS 
stars determined to be at that focus value. There 
is fair quantitative agreement between the data 
and the model. This agreement is poor in the 
center of the ACS WFC field. We have deter- 
mined that part of this disagreement is due to 
CTE degradation, which is not included in the 
TinyTim PSF models. We discuss the causes of 
this and our solution in § 15.21 the shown TinyTim 
model represents the best match to the ACS COS- 
MOS data from all the focus values for which we 
calculated TinyTim focus models. The real test 
of the quality of our PSF models is in the star- 
galaxy correlation functions before and after PSF 
correction. These are shown in Figure [TSl 

We have tested our focal position determina- 
tion using randomly selected stars taken from the 
calibration images of 47Tuc described in § 12.11 We 
select 20 random non-overlapping sets of 20 stars 
from the 746 stars seen in one of the 47Tuc images. 
We determine the focus for each of these sets of 
stars. As these stars are all from the same image, 
the focal position determined should be identical 
for all sets. We find that all of the determined fo- 
cal positions are within ±2/im of the mean focal 
position. Based on these simulations, the results 
of which are shown in Figure [9l we estimate the 
uncertainty with which we can measure the focal 
position in a short exposure to be about ±l/im. 
The uncertainty in determining the focus value for 
a COSMOS image is somewhat larger than this be- 
cause we average over the intra-orbit focus changes 
caused by telescope 'breathing' (see Figure [2]). 

5. PSF Correction in COSMOS Images 

We follow the procedure given in RRG (2000) 
to correct galaxy shapes for PSF. Stellar mo- 
ments are used to correct galaxy moments first 
for the isotropic portion of the PSF, then for the 
anisotropic portion. Quantities are kept in terms 
of moments (rather than ellipticities) through the 
entire correction process and only then are the cor- 
rected moments used to calculate galaxy elliptic- 
ities, and, ultimately, shear. We tested using a 
variable stellar weight function width w instead of 
the standard (for RRG) fixed stellar weight func- 
tion width and found that this did not significantly 
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Fig. 8. — A TinyTim PSF model (left panel) for 
a focus value of — 2/im and the average of many 
observed stars (right panel) from COSMOS fields 
with a similar apparent focus. There is rough 
qualitative agreement between the data and the 
models over much of the ACS field. The center 
of the chip does not show good agreement due, at 
least in part, to the effects of degradation of the 
CTE. The real data show less positive ei (elon- 
gation along the x axis) than the models, consis- 
tent with a CTE-induced smearing in the y di- 
rection. We remove this instrumental signature 
late in our lensing analysis as described in ^5.2\ 
Notice that the x and y ranges are not [0,4096] 
as in default ACS images but are larger because 
we are using a smaller output pixel scale. These 
plots are not meant to demonstrative quantitative 
agreement between the models and the real stars. 
The TinyTim models are noiseless, while the real 
COSMOS stars have image noise and each 'tick 
mark' represents a different number of stars (usu- 
ally about 10 or fewer). For the RRG method, 
it is important for the TinyTim stars to match 
the moments (rather than the more easily plot- 
ted ellipticities) of the real stars. A quantitative 
analysis of how well the model stars can be used 
to deconvolve the PSF can be obtained from the 
star-galaxy correlation functions in Figure [15] and 
in Massey et al. 2006a. 



improve the PSF correction. The weight function 
used to evaluate galaxies still varies with the size 
of the galaxy as discussed in § 12.41 

5.1. Applying the TinyTim models 

For each COSMOS field, we choose the best fit 
focus value for the telescope as described above. 
For each galaxy in that field we correct the mea- 
sured galaxy moments for PSF effects with the 
TinyTim model star at that focus which is clos- 
est to the galaxy's position. We have obviated the 
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Fig. 9. — Determination of the focal position from 
20 sets of 20 independent stars, each chosen ran- 
domly from a single image taken of a stellar field 
near the globular cluster 47Tuc. The x-axis shows 
the set number, from 1 to 20, and the y-axis shows 
the determined focal position. The uncertainty on 
the y-axis is the measured uncertainty of the fo- 
cal position. As is shown, the focal position de- 
termination is highly repeatable in random trials 
of independent stars at identical focal positions. 
From these sorts of tests, we estimate that the un- 
certainty in determining the focal position for a 
short exposure is approximately ±l/zm. The de- 
termination of the focus from a a set of exposures 
taken over a full orbit has a greater uncertainty 
due to the intra-orbit 'breathing' of the telescope. 
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Fig. 10. — The left panel shows the values of 
the Ixx moments for COSMOS stars and TinyTim 
model stars as a function of telescope focus val- 
ues, there is good agreement between the two. 
However, the right hand panel shows that the 
TinyTim models are not perfect. The fourth or- 
der moment Ixxxx is consistently underestimated 
by TinyTim. The correction factors between the 
real and TinyTim moments are tabulated in Ta- 
ble O 



as a perturbative correction to the second order 
moments in RRG, so this slight discrepancy be- 
tween model and real stars does not significantly 
effect our PSF correction. We verified this by mul- 
tiplying the TinyTim fourth order moments by a 
correction factor and re-running our PSF correc- 
tion routines and found no difference in residual 
PSF systematics. For completeness, we tabulate 
the multiplicative factor between the TinyTim and 
COSMOS stars in Table [51] It is unclear why the 
TinyTim models have fourth order moments that 
are slightly too small, but it may have to do with 
how the charge diffusion kernel is applied. 



need to interpolate the PSF across the field by cre- 
ating sufhciently dense grids of model PSFs. We 
have found that the TinyTim models match the el- 
lipticites (and more importantly the second order 
moments) of the COSMOS stars they are meant 
to represent. However, the TinyTim models do not 
have the exact profile of the real COSMOS stars 
and thus the fourth order moments are somewhat 
too small in the TinyTim models. Only the mo- 
ments with even powers of both x and y are af- 
fected. We show, for a range of focus values, the 
average of Ixx and the average of Ixxxx for both the 
TinyTim models and the COSMOS stars in Fig- 
ure 1101 The fourth order moments are used only 



5.2. Correction for Charge Transfer Effi- 
ciency (CTE) Degradation 

Gradual damage to CCD detectors due to expo- 
sure to charged particle radiation in the harsh en- 
vironment of space results in a degradation in the 
efficiency of charge transfer in the CCDs. High 
energy charged particle hits create charge traps 
that accumulate in the silicon substrate. These 
traps capture electrons for a short time, and re- 
lease them after a time delay drawn from an ex- 
ponential distribution. Several species of traps can 
exist, each with their own exponential time con- 
stant. When an image is read out after exposure, 
the electrons from sources are stepped across these 
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Ixxxx 


1.2 


^xxxy 
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Ixxyy 


1.1 


Ixyyy 


1 




1.2 



Table 1: Multiplicative factor between the size of 
moments in the TinyTim models and measured 
values of the COSMOS stars. This is the factor 
that the TinyTim models would have to be multi- 
plied by to have them equal in size to the moments 
of the COSMOS stars. 



charge traps. Their capture and subsequent re- 
lease trails the electrons away from the original 
sources, across the image in the opposite direc- 
tion to the readout. This spreading of the charge 
in the readout direction creates problems for pho- 
tometry, astrometry, and galaxy shape measure- 
ment. Significant effort has already gone into un- 
derstanding the CTE of the ACS WFC (Mutchler 
& Sirianni 2005) and correcting photometric errors 
due to charge loss (Riess & Mack 2004). However, 
the subtle changes in the shapes of objects due 
to CTE degradation have not yet been adequately 
addressed. 

Since ACS was installed in 2002, charge traps 
have accumulated to the level of tens per pixel. 
There are at least three different trap species, each 
with different release times for trapped electrons 
(Sirianni 2006). On the two-CCD WFC, charge 
readout occurs at the top of the top chip and the 
bottom of the bottom chip; charge is read out from 
the center of the field to the top and bottom edges. 
Charge incident furthest from the readout regis- 
ters is transferred over more pixels and therefore 
interacts with more traps during its translation 
to the chip edge. The effect on the shape of any 
one galaxy is difficult to predict, as it non- linearly 
depends upon the galaxys magnitude, size, ellip- 
ticity and radial profile, as well as its position on 
the CCD. Furthermore, the fixed number density 
of charge traps has the consequence of affecting 
faint sources more than bright ones; a hundred 
delayed electrons are significant in a source con- 
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Fig. 11. — An estimate of the spurious elliptic- 
ity induced by CTE charge trailing in a barely 
resolved, circular galaxy as a function of galaxy 
flux and at various positions on the CCD. In prac- 
tice, the actual amount of spurious ellipticity also 
depends upon the intrinsic ellipticity and the the 
radial profile of the object. 



taining a thousand electrons, but not one contain- 
ing 100,000. Thus, the shapes of distant galaxies 
are smeared by this effect, but less so the bright 
stars that are typically used to calibrate and test 
the galaxy shape measurement algorithms neces- 
sary for weak lensing. Thus, CTE degradation is 
a particularly difficult systematic effect to correct, 
because it cannot be calibrated using bright, high 
S/N stars as is typically done with other system- 
atics that affect weak lensing measurements. 

The effect of the CTE degradation is a smear- 
ing of objects in the readout (y) direction. For 
faint galaxies, the size of this effect is comparable 
to the size of the weak lensing signal we are trying 
to detect, as shown in Figure [TTl The smearing is 
caused by each pixel leaving an exponential trail 
during readout as shown in Figure [T2l We have de- 
veloped preliminary models for what CTE degra- 
dation does to the charge in individual pixels and 
we use that to show the effect of CTE degradation 
on the image of a faint galaxy in Figure 1131 

Previous weak lensing work with HST has en- 
countered similar problems. Rhodes et al. (2004) 
corrected for CTE in STIS by creating S/N de- 
pendent PSF models. That is impractical here be- 
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Fig. 12. — Exponential profile created during 
charge readout as measured in an EPER image 
as described in § 12.11 The value of the pixel that 
this tail corresponds to is about 2500 electrons. 



cause we are already dealing with focus-dependent 
PSF models and the computing power and com- 
plexity required to model stars across the ACS 
WFC field at a range of focus values and S/N 
values would quickly become prohibitive. Further- 
more, the CTE degradation grows worse with time 
and the COSMOS images are taken over a rela- 
tively long time period. Ideally, we would correct 
each pixel of each image for the effects of CTE 
degradation as the first step in the image reduc- 
tion pipeline (see for example Bristow et al.2002). 
CTE effects are the last to go into the image since 
they occur during the readout. However, this 
would require an exact knowledge of the number of 
charge traps, the number of species of charge traps 
(there are at least 3 species each with different re- 
lease times), and accurate knowledge of the release 
times of those charge traps. We are working on a 
general solution to this problem that incorporates 
all of that knowledge. In the meantime, we have 
developed a parametric equation that allows us to 
correct the measured ellipticities of galaxies for the 
effects of CTE degradation. The CTE depends on 
the position of the object within the CCD (elec- 
trons further from the readout registers encounter 
more charge traps), the flux of an object (high flux 
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Fig. 13. — Illustration of the effect of charge 
trailing on a 25th magnitude star in single-orbit 
F814W data, using a crude pixel- level CTE model 
we have developed. The readout direction is down. 
A simulated galaxy is shown in the left panel as 
it ought to appear and, in the right panel, after 
being read out from the furthest side of the CCD. 
The difference between the images is small and 
cannot be detected by eye in the images shown. 
The difference is, however, significant at the level 
of precision needed for weak lensing. This differ- 
ence is shown in the bottom panel (with a greatly 
expanded grey scale). We eventually hope to un- 
derstand how charge is transferred between indi- 
vidual pixels in the ACS CCD. For this work we 
correct for the overall effect the CTE has on galaxy 
ellipticities with Equation fTUl 

objects fill the charge traps and the relative loss of 
flux is less), and the date of observation (the CTE 
is continually degrading due to cosmic ray dam- 
age). By assuming that the PSf corrected elliptic- 
ities of all the galaxies in the COSMOS fields aver- 
age to zero, Massey et al. 2006a have found empiri- 
cally the dependence of CTE effects on these three 
variables and we use this empirical knowledge to 
derive the parametric correction equation 
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= ef - eete (9) 

where el is the CTE corrected first elhpticity com- 
ponent in the undistorted image coordinate sys- 
tem (which is changed by only a few degrees to 
the distorted coordinates), and ef is the mea- 
sured first eUipticity component after correction 
for other PSF affects using the RRG method, as 
described above. Only the first elhpticity compo- 
nent is affected because this component represents 
elongation in the x and y direction. 62 represents 
elongation along axes at ±45° to the x axis. The 
additive correction parameter e^^g is given by 

e^te = -3xl0-5xS/N"\^^^^g^^x(MJ£i-52333) 

(10) 

where S/N is the object's detection signal-to-noise, 
ntransfors is the distancc to the nearest readout reg- 
ister in native ACS pixels, and MJD is the modi- 
fied Julian date of the observation. Note that this 
empirically-derived correction depends on date, y 
position (number of readouts), and fiux, just as 
the CTE effect is known to. We were able to ob- 
tain similar results by adjusting the exponent on 
the S /N term while adding a size-dependent term. 
However, we choose this formalism because it is 
simpler and more physically motivated. It is im- 
portant to note that this particular formalism for 
the CTE is only valid on data taken before July 
2006. At this time the operating temperature of 
the ACS WFC CCDs was changed, thus chang- 
ing the CTE of the CCDs (Sirianni, GiUiand, & 
Sembach 2006). 

We show in Figure [14] the ellipticities of galax- 
ies as a function y position and magnitude both 
before and after the CTE correction for the en- 
tire COSMOS survey The tell-tale dip in ef , 
implying an elongation in the y direction, found 
in the data is removed by this parametric CTE 
correction. Before correction, the mean elhptic- 
ity of the whole galaxy population is {ef') = 
-0.020 ± 0.001, and the gradient of a fit to the 
faintest magnitude bin in Figure [HI is „ ^'^^ — = 

CTltransfcra 

(-1.74 ± 1.3) X lO"'^. After correction, these 
values are reduced to (e^) = 0.004 ± 0.002 and 

^5-^£J = (-1.53 ± 1.0) X lO-'^. Note that in 

two-point shear correlation functions (e.g. Massey 
et al. 2006a) this residual enters only as the value 



squared. In the sense of that statistic, we have 
therefore lessened the impact of CTE trailing in 
faint galaxies by more than two orders of magni- 
tude. 

Despite its apparent success, we stress that this 
prescription is by no means a panacea for CTE 
effects in ACS weak lensing data. In this sim- 
plified model, we take advantage of the uniform 
background level of the COSMOS images to elim- 
inate dependence on this parameter. The other 
model parameters are also specific to our dataset 
and shear measurement method. An improved, 
pixel-level CTE correction method, along the lines 
of Bristow et al. (2002) will follow (Massey & 
Rhodes, 2007). This method will take into account 
the different species of charge traps and their as- 
sociated release times. 

5.3. Performance of the PSF Correction 

To demonstrate the correction of the ACS data 
for PSF and CTE degradation, we show correla- 
tion functions between the galaxy shears and the 
raw stellar ellipticities in Figure [151 Stars and 
galaxies should be correlated before PSF correc- 
tion (due to convolution by the same PSF) but 
should be uncorrelated after PSF correction, note 
that in the RRG method we do not correct the 
stars for PSF convolution. Thus, after PSF cor- 
rection the stars should contain the PSF signal 
but the galaxies should not, and their ellipticities 
should be uncorrelated. Further descriptions of 
the use of correlation functions in weak lensing 
can be found in, for instance. Bacon et al.2003 and 
Kamionkowski et al.l998. This figure shows 

cm = ^{erir)Yiir + 0)) , (11) 

(Te* 

C2iO) = ^ ( er(r) 72^(r + 0) ) and (12) 

(7 e* 

Gq* I 

(er(r) lliv + O))^ ,(13) 

where 6 is the separation between a star and the 
galaxies, averaging is performed over the whole 
population and the superscript ^ denotes compo- 
nents of elhpticity, and shear rotated so that 7]" 
(72) ill each galaxy points along (at 45° from) the 
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Fig. 14. — The top panel shows the average elhp- 
ticity as a function of y position for several magni- 
tude ranges after the RRG PSF correction scheme. 
The lowest line in the top plot represents galaxies 
in the magnitude range 26-27. The other lines rep- 
resent one magnitude bins going up to 22-23. The 
effects of CTE are clearly shown by the fact that 
fainter objects (with lower flux) are preferentially 
elongated in the y direction, corresponding to neg- 
ative ei- The problem is worst near the center of 
the field, farthest from the readout registers. The 
bottom panel shows that a simple additive fac- 
tor (as given by Equation [TOl) on the ellipticity of 
each galaxy removes the CTE-induced ellipticity. 
In the future we hope to understand the CTE well 
enough to remove CTE effects from the images 
before a lensing analysis. 

vector between the pair. The normalization via 
star-star correlation functions suggested by Bacon 
et al. (2003) is impractical as a denominator in 
this case, because the specific PSF pattern of ACS 
makes it cross zero several times. Figure [13] also 
shows the star-galaxy correlation functions after 
separation into E-l-B modes via the variance of 
the aperture mass statistic, as defined in Schnei- 



der, van Waerbeke & Mellier (2002). This figure 
includes error bars from statistics alone (inner) 
and, when applicable, variation found be subdi- 
viding the COSMOS field into four independent 
quadrants and doing the analysis on each quad- 
rant independently (outer). These error bars in- 
crease in size after the PSF correction. This is 
due to the fact that PSF convolution circularizes 
objects (making their ellipticities smaller). After 
PSF correction, the ellipticity distribution of ob- 
jects is larger and thus the scatter in the correla- 
tion functions is larger. 

Before correction, the shear measurements con- 
tain artifacts from the PSF anisotropy. These are 
largely removed by the processes described in this 
section, and are consistent with zero on all scales 
after correction. Correlation function Ci and C2 
show significant improvement (movement towards 
zero correlation) , especially on scales less than one 
arcminute. C3 was already nearly consistent with 
zero before PSF correction and the PSF correction 
has introduced a larger scatter due to the widen- 
ing of the ellipticity distribution during PSF de- 
convolution discussed in the previous paragraph. 
The most important test of our PSF correction is 
shown in the E-l-B mode plots, which show that 
after PSF correction the E-fB mode is consistent 
with zero for all scales. See Massey et al. 2006a 
for a more detailed description of these correla- 
tion functions and the separation of the signal into 
E+B modes. 

6. Conclusions 

We have shown that aliasing of the ACS WFC 
PSF can be minimized by carefully choosing the 
parameters of the image reduction pipeline ele- 
ment MultiDrizzle. This aliasing is minimal 
when MultiDrizzle is run with a Gaussian ker- 
nel, the area of the input pixels containing the 
flux is shrunk by a factor of pixfrac= 0.8, and 
the output image has a pixel scale of 0.03 arc- 
seconds per pixel. We show that the ACS WFC 
PSF is temporally unstable over the time scale of 
the COSMOS observations due to thermally in- 
duced changes in the telescope's focus. Using a 
modified version of the TinyTim software package 
we create dense PSF grids at a range of telescope 
focus values from — 10/xm to -|-5/im. Using the 
~ 10 suitable stars in each image as taken from the 
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Fig. 15. — Correlation functions between galaxy 
shears and raw sellar ellipticities, showing the con- 
tamination of the shear catalog by residual PSF 
artifacts. The inner error bars contain statisti- 
cal errors only and the outer error bars include 
the variation found by dividing the COSMOS field 
into four quadrants and performing the correlation 
function analyses on each quadrant independently. 
The left hand panels show the shear measured 
straight from the images; the right hand panels 
show the shear after correction. The solid lines 
show the best-fitting cosmic shear signal measured 
by Massey et al. (2006a) for comparison. The suc- 
cess of our correction scheme is demonstrated by 
the lowering of the data points to be consistent 
with zero after correction, especially in the plot of 
E and B mode signal. 

COSMOS lensing catalog (Leauthaud et al. 2006), 
we can calculate the focus value of HST for each 
COSMOS field within about ±l/xm. We can use 
the TinyTim models to correct the galaxy shapes 
for the effects of PSF. There is a residual system- 
atic left after the PSF correction that is due to 
degradation of the CTE by cosmic ray damage to 
the CCDs. We present a parametric correction 
for the effects of CTE degradation on galaxy el- 



lipticities. The resulting PSF and CTE corrected 
catalog has been used to measure cosmic shear by 
Massey et al. (2006). 

Our wrappers and implementation of the 
TinyTim code are publicly available from 
\protectlittp://www. astro. caltech.edu/~rjm/acs 
Since our PSF modelling routines are general, they 
are applicable to other filter sets and can be used 
for weak lensing data sets taken with other fil- 
ters (e.g. the ACS parallel survey taken with the 
F775W filter or a re-reduction of the GEMS sur- 
vey which was taken with the F606W filter). Sim- 
ilarly, the PSF models are not specific to the RRG 
weak lensing method. Since we are creating sim- 
ulated stars, any weak lensing pipeline can make 
use of these stars for PSF deconvolution. The code 
and models have already been put to use in other 
science papers by the COSMOS collaboration in 
such diverse areas as AGN studies (Gabor et al. 
2006) and morphological classification (Scarlata 
et al. 2006). 

Future work will include a more general cor- 
rection for the effects of CTE degradation. This 
correction will likely take place on the images 
themselves as a first step in the image reduction 
pipeline. Given the general nature of our PSF 
models, we plan to use them to perform further 
weak lensing analyses on the COSMOS data set 
with the more advanced "shapelets" method (Re- 
fregier 2003; Refregier & Bacon 2003; Massey & 
Refregier 2005). 
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